Datos sobre primates.

inp<- 
  read.csv("primates-cr-registros.csv", sep = "\t")

Datos de registros en la Tabla.

inp %>%
  dplyr::select(family,
         species,stateProvince,
         locality,eventDate) %>%
  DT::datatable(
    colnames = c("Familia", "Especie","Provincia","Localidad", "Fecha"),
    rownames = FALSE,
    options = list(
      searchHighlight = TRUE,
      language= list(url = '//cdn.datatables.net/plug-ins/1.10.11/i18n/Spanish.json'),
      pageLength = 15))

Grafico de pastel sobre registros de las Especies

primates <- inp[,c('species', 'eventDate')]
data <-primates %>%
  group_by(especies = format(species)) %>%
  summarize(suma_registros = n())

data$suma_registros<- as.numeric(data$suma_registros)

plot_ly(data, labels = ~especies, values = ~suma_registros, 
        type = 'pie') %>% config(locale= "es") %>% 
  layout(title = "Grafico Registros de Especies")

#Filtro de datos

AP<-filter(inp, species == "Alouatta palliata")
CC<-filter(inp, species == "Cebus capucinus")
AG<-filter(inp, species == "Ateles geoffroyi")
S<-filter(inp, species == "Saimiri oerstedii")

MAPA

#Lectura de una archivo CSV con columnas de coordenadas
inp1 <-
  st_read(
    "primates-cr-registros.csv",
    options = c(
      "X_POSSIBLE_NAMES=decimalLongitude",
      "Y_POSSIBLE_NAMES=decimalLatitude"),
    quiet = TRUE)

# Asignación del sistema de coordenadas

st_crs(inp1) = 4326
# Mapeo
mapa <- inp1 %>%
  dplyr::select(family,
         species, stateProvince,
         locality,
         eventDate) %>%
  leaflet() %>%
  addProviderTiles(providers$OpenStreetMap.Mapnik, group = "OpenStreetMap") %>%
  addProviderTiles(providers$Esri.WorldImagery, group = "Imágenes de ESRI") %>%
  addCircleMarkers(
    stroke = F,
    radius = 2,
fillColor = 'red',
    fillOpacity = 1,
    popup = paste(
      inp$family,
      inp1$species,
      inp1$stateProvince,
      inp1$locality,
      inp1$eventDate,
      sep = '<br/>'),
    group = "species") %>% addLayersControl(
      baseGroups = c("OpenStreetMap", "Imágenes de ESRI"),
      overlayGroups = c("primates")) %>%
  addMiniMap(
    tiles = providers$Stamen.OpenStreetMap.Mapnik,
    position = "bottomleft",
    toggleDisplay = TRUE)
# 
alt <- getData(
    "worldclim",
    var = "alt",
    res = .5,
    lon = -84,
    lat = 10)
  
  altitud <- crop(alt, extent(-86, -82.3, 8, 11.3))
  
pal <- colorNumeric(
  c("#006400", "#FFFF00", "#0000FF"), 
  values(altitud), 
  na.color = "transparent")
mapa %>%
  addRasterImage(
    altitud, 
    colors = pal, 
    opacity = 0.8) %>% addLayersControl(
      baseGroups = c("OpenStreetMap","Imágenes de ESRI","Worldclim"))
## Warning in showSRID(uprojargs, format = "PROJ", multiline = "NO", prefer_proj =
## prefer_proj): Discarded ellps WGS 84 in Proj4 definition: +proj=merc +a=6378137
## +b=6378137 +lat_ts=0 +lon_0=0 +x_0=0 +y_0=0 +k=1 +units=m +nadgrids=@null
## +wktext +no_defs +type=crs
## Warning in showSRID(uprojargs, format = "PROJ", multiline = "NO", prefer_proj =
## prefer_proj): Discarded datum World Geodetic System 1984 in Proj4 definition
## Warning in showSRID(uprojargs, format = "PROJ", multiline = "NO", prefer_proj =
## prefer_proj): Discarded ellps WGS 84 in Proj4 definition: +proj=merc +a=6378137
## +b=6378137 +lat_ts=0 +lon_0=0 +x_0=0 +y_0=0 +k=1 +units=m +nadgrids=@null
## +wktext +no_defs +type=crs
## Warning in showSRID(uprojargs, format = "PROJ", multiline = "NO", prefer_proj =
## prefer_proj): Discarded datum World Geodetic System 1984 in Proj4 definition